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Abstract. Schlesinger transformations are algebraic transformations of a Fuchsian system that preserve its 
monodromy representation and act on the characteristic indices of the system by integral shifts. One of the 
important reasons to study such transformations is the relationship between Schlesinger transformations and 
discrete Painleve equations; this is also the main theme behind our work. In this paper we show how to write 
an elementary Schlesinger transformation as a discrete Hamiltonian system w.r.t. the standard symplectic 
structure on the space of Fuchsian systems. We then show how Schlesinger transformations reduce to 
discrete Painleve equations by considering two explicit examples, d-P(Z)^^') (or difference Painleve V) and 
d-P(AW*). 



1. Introduction 

In the theory of ordinary hnear differential equations on a complex domain, and in particular in the 
theory of Fuchsian systems, an important characteristic of the equation is its configuration of singularities 
and the characteristic indices at these singular points. Associated to this data is the notion of the monodromy 
representation of the equation. Roughly speaking, this representation describes how the fundamental solution 
matrix of the equation changes under analytic continuation along the closed paths around the singular points 
and it gives a significant insight into the global behavior of solutions of the equation. 

The fruitful idea of deforming the equation by moving the location of critical points into different con- 
figurations, or by changing the characteristic indices, without changing its monodromy representation, goes 
back to B. Riemann, but the actual foundations of this theory of isomonodromic deformations in the Fuch- 
sian case were laid down in the works of R. Fuchs [Fuc07], L. Schlesinger [Schl2], and R. Garnier [Gar26]. 
An extension of the theory to the non-Fuchsian case was done relatively recently in the series of papers by 
M. Jimbo, T. Miwa, and K. Ueno, [JMU81, JM82, JM81]. At present, the theory of isomonodromic defor- 
mations is a very active research field with deep connections to other areas such as the theory of integrable 
systems, classical theory of differential equations, and differential and algebraic geometry. 

In this work we focus on the relationship between isomonodromic deformations and Painleve equations. 
We need to distinguish between continuous and discrete isomonodromic deformations. In the continuous 
case we move the location of singular points and the resulting Schlesinger equations reduce to Painleve-type 
nonlinear differential equations. On the other hand, if we deform the characteristic indices, the isomon- 
odromy condition requires that the indices change by integral shifts, and so the resulting dynamic is discrete 
and is expressed in the form of difference equations called Schlesinger transformations. Similarly to the con- 
tinuous case, these transformations reduce to discrete analogues of Painleve-type equations and the study 
of this correspondence has been a major research topic in the field of discrete integrable systems over the 
last twenty years, [RGII91, PNGR92]. We need to remark here that Schlesinger transformations correspond 
to difference Painleve equations, but the other two types of discrete Painleve equations, q-difference and 
elliptic- difference Painleve equations can also be considered in a modification of this approach. 

Discrete Painleve equations share many properties with the differential Painleve equations, e.g., the 
existence of special solutions such as algebraic solutions or solutions that can be expressed in terms of special 
functions, affine Weyl group symmetries, and the geometric classification of equations in terms of rational 
surfaces. However, while it is well-known how to write differential Painleve equations in the Hamiltonian 
form [OkaSO], the discrete Hamiltonian structure for discrete Painleve equations is at present missing. 
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The aim of the present paper is to study such discrete Hamiltonian structure for Schlesinger transforma- 
tions. In particular, we present an exphcit formula for the discrete Hamiltonian function of an elementary 
Schlesinger transformation expressed in the same coordinates as the Hamiltonian function of the continuous 
Schlesinger equations considered by Jimbo, Miwa, Mori, and Sato in [JMMS80]. Note that this explicit for- 
mula gives a convenient tool for computing Schlesinger transformations and consequently for the derivation 
of discrete Painleve equations, since the usual computation of a Schlesinger transformation from the com- 
patibility condition between the Fuchsian equation and the deformation equation is often quite complicated. 

We also present two explicit examples of discrete Painleve equation computed in this framework. The first 
example that we consider is an isomonodromic deformation of a two-dimensional Fuchsian system with three 
finite singular points. A nice feature of this example is that it allows us to compare, in the same setting and 
using the same coordinates, continuous deformations that are described by the Painleve- VI equation Pyi, 
and discrete deformations that are described by the difference Painleve- V equation d-Py, which corresponds 
to the Backlund transformations of Pyi- Recall that in terms of the surface of initial conditions d-Py is also 
denoted as d-P{D';^''). In the second example we consider a discrete Painleve equation of the type d-P(v42^''*) 
that does not have a continuous Painleve counterpart. In [Sak07] one of the authors posed a problem of 
representing this equation using Schlesinger transformations of some Fuchsian system. In [Boa09] P. Boalch 
described the Fuchsian system whose Schlesinger transformations are described by d-P(^2^''*) without doing 
the explicit coordinate computation. Also note that D. Arinkin and A. Borodin clarified the correspondence 
between d-P^A^'^*) and the difference Fuchsian equation in [AB06]. Here we compute the d-P(742^-'*) 
example explicitly using the discrete Hamiltonian and use it to illustrate some interesting features, as well 
as some difficulties, that occur when we study discrete Painleve equations using this approach. 

The text is organized as follows. In the next section we briefly set up the context of the problem and 
describe in more details some of the main objects. In Section 3 we derive the equations of the elementary 
Schlesinger transformations and in Section 4 we show how to write them in the Hamiltonian form. Section 5 
is dedicated to the d-P{D'^^) and d-P(^2^''*) examples, and the final section is conclusions and discussion 
of the resul. 
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2. Preliminaries 

2.1. Fuchsian Equations. Schlesinger transformations that we study originate from the deformation theory 
of Fuchsian equations. Recall that a Fuchsian equation (or a Fuchsian system) is a matrix linear differential 
equation on the Riemann sphere CP such that all of its singular points are regular singular points. We 
consider a generic case when a Fuchsian equation can be written in the the Schlesinger normal form, i.e., 
when the coefficient matrix is a rational function with at most simple poles at some (distinct) points wi, . . . , u„ 
(and possibly at the point uq — oo), 

Here = res„. A.[x) dx are constant m x m matrices. We also put Aqq — le&oa A(a;) dx — ~ ^'^'^ 
so this equation is regular at infinity iff Aqo = 0. 

Assumption 2.1. From now on we make an additional semi-simplicity assumption that the coefficient 
matrices A^ are diagonalizable. 

2.2. Spectral Type and Accessory Parameters. Geometrically, Schlesinger (or isomonodromic) dy- 
namic takes place on the space of coefficients A{x) of the Fuchsian equation (2.1), considered modulo gauge 
transformations. Here we briefly outline the description of this space following [SaklO]. 
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First, given A(a;), we separate the data about locations of singular points Ui (that we think of as parameters 
of the dynamic) and the residue matrices at those points. Thus, we define 

^(u) = ^(ui, . . . ,M„) - {(Ai, . . . , A„) I A, e End(C™)} 

to be the set of all Fuchsian equations with possible singularities at the points Ui. Further, since Schlesinger 
dynamic either preserves (in the continuous case) or shifts (in the discrete case) the eigenvalues of A^, we 
treat the eigenvalues as the parameters of the dynamic as well. Thus, the appropriate phase space is a 
quotient of the fiber of the eigenvalue map by gauge transformations. Local coordinates on the phase space 
are called accessory parameters. Note that the space A of accessory parameters is quite complicated and its 
dimension depends on the spectral type of A(a;) that we define next. 

Spectral type of A(a::) encodes the degeneracy of eigenvalues of the coefficient matrices A^ via partitions 
m = m\ + • • • + m[% m\ > • ■ ■ > > 1, where m is the matrix size and denotes the multiplicities of 
the eigenvalues of A^. The spectral type m of A(x) is then defined to be the collection of these partitions for 
all indices (including oo): 

m = • • • m'lS m2 • • • • • • , TO^ • • • mj,", • • • . 

It can be used to classify Fuchsian systems up to isomorphisms and the operations of addition and middle 
convolution introduced by N. Katz [Kat96], as in the recent work by T. Oshima [Osh08]. 

Assumption 2.2. Since for a Fuchsian system we can use scalar local gauge transformation of the form 
Y(x) — z{x)^^\{x), where z{x) is a solution of the scalar equation 

dz el 

dx ^ X — Ui ' 

2—1 

to change the residue matrices by A^ = A; — 6*^1, we can assume, without loss of generality, that one of the 
eigenvalues d\ = 0. Thus, we always assume that the eigenvalue 6\ of the highest multiplicity m\ is zero 
and we denote the corresponding subset of T{\\) of reduced Fuchsian equations by -^(u). This amounts to 
choosing a representative in the quotient space of all Fuchsian equations by the group of local scalar gauge 
transformations. 

In view of Assumptions 2.1 and 2.2, A^ is similar to a diagonal matrix diag{6'|, . . . , 6'[' , 0, . . . , 0}, where 
ri = rank(Aj). Omitting the zero eigenvalues, we put 

0, = diag{0i,...,0r}. 

Denote by T the set of all possible diagonal matrices of the spectral type m with the highest degeneracy 
eigenvalues at finite points set to and omitted, as above. Thus, T = {0 = (0i, . . . , 0„; ©oo)} and 
ri = m — mj. Then we have an eigenvalue map Sp : J-{u) — >■ T from the set of reduced Fuchsian equations 
to the set of characteristic indices, where T C T is the subset of © satisfying the Fuchs relation (or the 
trace condition) 

n,oo n 

i=i j=i 

With these definitions, given © S T, the fiber of the eigenvalue map is J^(u)© = Sp~"'^(©). 

We still need to take into account the global similarity transformations. It is convenient to do this in 
two steps. First, we use such transformations to normalize our equation at infinity by reducing Aqo to a 
particular form, e.g., make it diagonal, and then take the further quotient by the stabilizer subgroup Ga^o 
of the group GL™ of global similarity transformations. Thus, for a given © G T, we fix Aoo such that 
Sp(Aoo) = ©oo and denote by J^(u)(0_Aoo) the subset of all Fuchsian equations in ^^(u)® satisfying the 
condition Ai + • • • + A„ = — Aoo- Then, finally, our phase space is 
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When Aoo is diagonal, which is the case that we are mostly interested in, the dimension of this space A 
of accessory parameters of the spectral type m is given by the formula 




dimyl= (n- 1)to^ - > >(m^') +2. 



2.3. The Decomposition Space. To study the Hamiltonian structure of Schlesinger equations, and also 
of Schlesinger transformations, it is convenient to consider a larger decomposition space that is defined as 
follows. 

In view of Assumption 2.1, there exist full sets of right eigenvectors hij, Aihij — O-fhij, and left 
eigenvectors cj\ cj^ -^i = ^icf' , (here we use the f symbol to indicate a row-vector). In the matrix form, 
omitting vectors with indices j > ri that are in the kernel of A^, we can write 



B, = [b,,i • • • b,,,J , A,B, = B,0„ C| - 



C A,; - 0,C 



with &i defined by (2.2). But this means that we have a decomposition A; — B^Cj provided that CjBi = 0;, 
where the last condition is related to the normalization ambiguity of the eigenvectors. 

Thus, given A^, we can construct, in a non-unique fashion, a corresponding decomposition pair (B.^, C|^). 
The space of all such pairs for all finite indices 1 < i < n, without any additional conditions, we call the 
decomposition space. We denote it as 

BxC = (C' X • • • X C") X {{C^y X • • • X (C")^) ^ (C^ X (C^Y) X • • • X (C^" x (C''")^), 

since it is convenient to write an element (B, C^^) of this space as a list of n pairs (Bi, C|; • • • ; B„, CjJ. This 
space can be equipped with the natural symplectic structure where we take the matrix elements of B^ and 
as canonical coordinates, i.e., we take our symplectic form to be 

n 

a; = ^tr(dCj AdB,). 

1=1 

Remark 2.3. To avoid symbol overcrowding we slightly abuse the notation and use the same symbols, e.g., 
Bi or bi.j, to denote both the canonical coordinate system in the decomposition space and an actual point 
in the space; the exact meaning of the symbol is always clear from the context. 

The idea of using the decomposition space B x C to study the Hamiltonian structure of isomonodromic 
deformations goes back to the paper by M. Jimbo, T. Miwa, Y. Mori, and M. Sato [JMMS80]. Since then it 
has been used by many other researchers studying isomonodromic deformations, most notably by J. Harnad, 
see, e.g., [IIar94]. 

Remark 2.4. There are two natural actions on the decomposition space B x C. First, the group GL„i of 
global gauge transformations of the Fuchsian system induces the following action. Given P G <GLm, we have 
the action A^ PA^P"^ which translates into the action (Bi,Cj) H> (PBi,CjP~^). We refer to such 
transformations as similarity transformations. It is often necessary to restrict this action to the subgroup 
Ga„ preserving the form of Aqo- Second, for any pair (Bi,Cj) the pair (BiQi,Q~^C^) determines the 
same matrix A^ for Qi € GL^^. The condition Q,^^C|BiQi = Q~^0iQi = &i restricts Qi to the stabilizer 
subgroup Ge; of Ghr-. In particular, when all are distinct, Qi has to be a diagonal matrix. We refer to 
such transformations as trivial transformations. These two actions obviously commute with each other. 

We are now ready to give an alternative description of the space of accessory parameters ^(q.Aoo)- Given 
the pair (0, Aoo) as in the previous section, let 

n 

{B X C)(0,A„) = {(Bi, Cl ; • • • ; B„, Cl) e B x C \ c|b, = 0„ ^ B,c| = -Aoo} cBxC. 



Then 



A = y^(e,A^) = X C)(©,A^)/((G'©i X • • • X G©J x Ga^). 
The following diagram illustrates the relationship between different spaces defined in these two sections: 

{B X C)(e,A^,) C B^iC 




-^(u)(e,A^) C J^(u) 

•A = ^(e,Aoo) 

Figure f . The decomposition space and the space of accessory parameters 



2.4. Schlesinger Equations and Schlesinger Transformations. Consider now isomonodromic defor- 
mations of a Fuchsian system. For the continuous case, the most interesting situation is when we take the 
deformation parameters to be locations of critical points. Thus, we let A = A(x, u) (i.e., A^ = Ai(u)), and 
consider a family of Fuchsian equations on the function Y = Y(x; u). 



A(x;u)Y= [J2 



dx \ ^ X — Ui 

In [Schl2] L. Schlesinger showed that the monodromy preserving condition is equivalent to the set of defor- 
mation equations 

dY A,(u) . 

— = Y, z = l,...,n. 

oui X ~ Ui 

The compatibility conditions for these deformation equations are called Schlesinger equations. These are 
partial differential equations on the coefficient matrices Ai(u) and they have the form 

OA, _ [Aj, Aj] _^ aA, _ _ [Aj, Aj] 

duj Ui — Uj ' dui ^-^ Ui — Uj 

In [JMMS80] Jimbo et.al showed that Schlesinger equations can be written as a Hamiltonian system on 
the decomposition space B x C, 



with the Hamiltonian 



duj 5ct ' duj dB 



n, ^ n,{B, ct) ^ A, = B.cl. 

Ui — Ui 



Schlesinger equations admit reductions to Painleve-type nonlinear differential equations, which remains valid 
on the level of the Hamiltonians as well; we review an example of Pyi in Section 5. The main question that 
we consider is what happens to this picture in the discrete case. 

The discrete counterparts of Schlesinger equations are Schlesinger transformations. Those are rational 
transformations preserving the singularity structure and the monodromy data of the Fuchsian system (2.1), 
except for the integral shifts in the characteristic indices Of . The study of such transformations again goes 
back to Schlesinger [Schl2], and this and the more general case of irregular singular points was considered in 
great detail in [JM81]. Schlesinger transformations have the form Y(a;) i-> Y{x) — R(a;)Y(2;), where 'R{x) 
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is a specially chosen rational matrix function called the multiplier of the transformation. The coefficient 
matrix A(a;) of om' Fuchsian system then transforms to A(x) that is related to A(x) by the equation 

dR{x) 



(2.2) 



A(a;)R(a;) = R(a;)A(x) 



dx 



Similarly to the continuous case, Schlesinger transformations induce discrete Painleve-type dynamic on 
the space of accessory parameters. This dynamic is very complicated, but it becomes much simpler when 
considered on the decomposition space, which allows us to represent it as a discrete Hamiltonian system 
using the same canonical coordinates as [ JMMS80] , as in Figure 2: 



BxC * BxC 

Discrete 



u 



(^xC)(e,A. 



A 



(®,A, 



Hamiltonian 
Dynamic 



Painleve 
Dynamic 



u 



^(©,a; 



Schlesinger 
transformations _ 

► 

r 



Figure 2. Schlesinger Transformations 



Unfortunately, at present we do not have a discrete version of the Hamiltonian reduction procedure from 
the decomposition space to the space of accessory parameters, this is an important problem for future 
research. 

In this paper we focus on the elementary Schlesinger transformations that only change two of the charac- 
teristic indices by unit shifts, i.e., 0^ = 0^ — 1 and 9^ =0^ + 1. For the examples section, we also assume 6*^ 
and 0^ have no multiplicities. In [JM81] such elementary transformation are denoted by { " }. In Section 4 
we construct a function 'H+(B, C^^) of the decomposition space BxC that is a discrete Hamiltonian function, 
in the sense explained in the next section, for such elementary Schlesinger transformation. 

2.5. Discrete Lagrangian and Hamiltonian Formalism. In developing a discrete version of the Hamil- 
tonian formalism, one approach is to discretize the continuous dynamic preserving its integrability properties, 
see the recent encyclopedic book by Y. Suris, [Sur03], but it requires underlying continuous dynamical sys- 
tem. An alternative and more direct procedure is to develop a discrete version based on the variational 
principles. In this case, the Lagrangian formalism is more natural, but it is possible to extend it to include 
the discrete version of the Hamiltonian formalism as well. This approach has its origin in the optimal control 
theory (see, e.g., B. Jordan and E. Polak [JP64], and J. Cadzow [Cad70]) and mechanics (see e.g., J. Logan 
[Log73] and S. Maeda [Mae82]). In the theory of integrable systems it was first used in the foundational 
works by A. Veselov [Ves91, VesSS] and A. Veselov and J. Moser [MV91]. A lot of recent work in this field 
is motivated by the development of very effective symplectic integrators, see an excellent review paper by 
J. Marsden and M. West [MWOl] (and the references therein), as well as more recent works by S. Lall and 
M. West [LW06] and A. Bloch, M. Leok, and T. Oshawa [OBLll] that emphasize the Hamiltonian aspects 
of the theory. Below we give a very brief outline of this approach closely following [OBLll], we refer to this 
and other publications above for details. 

Let Q be our configuration space. Then, in the discrete case, the Lagrangian C G J'iQ x Q) is a function 
on the state space Qx Q. For the discrete time parameter A: g Z, a trajectory of motion is a map q : Z — > Q, 
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q(fc) —: qfc, or, equivalently, a sequence {q^ S Q}. The action functional on the space of such sequences 
is given by ©({q^}) = £(qfc, qfc+i). Using the variational principle S& = we obtain discrete Euler- 
Lagrange equations that have the form 

L'2>C(qfc_i,qA;) + £'i£(qfc,qfc+i) = 0, 

where Di (resp. D2) denote vectors of partial derivatives w.r.t. the first (resp. second) sets of local coordinates 
on the state space Q x Q. These equations then implicitly determine the map (or more precisely, the 
correspondence) qt+i = 0(qfc_i, q^), which then defines the discrete Lagrangian flow FciQxQ— >QxQ 
on the state space via ^£(qfe-i,q*:) = (q/c^q/c+i)- This flow is symplectic w.r.t. the discrete Lagrangian 
symplectic form 

f^£(qfe,qfe+i) := dd]^ = d-d^ = L»ii:)2£(qfe, qfe+i)dqfc A dqfe+i, 
F^Q.c — ^c, where one-forms i?^ : Q x Q T* (Q x Q) are defined by 

(q/c, qfe+i) = D2£{qk, qfc+i)dqfc+i, ^9j(qfc, q/c+i) = -Di£{qk, qfc+i)rfqfc. 

If we introduce right and left discrete Legendre transforms ¥C^ : Q x Q ^ T* Q and the momenta 
variables pfc, p^+i by 

F/:+(qfc,qfc+i) = (qfc+i,Pfe+i) = (qfc+i, i:'2/:(qfc, qfc+i)), 
F/:"(qfe,qfe+i) = (qfc, Pfc) = (qfc, -L»i/:(qfc, qfc+i)), 

we see that -d^ — (F£*)*?? and i7£ — (F£*)*f2, where ?? and are the standard Liouville and symplectic 
forms on T* Q respectively. 

We can then define the discrete Hamiltonian flow Fn : T* Q T*Q by Fc(qfc,Pfc) = (qfc+i,Pfc+i), i-e., 
Fc = VC+ o (F£-)-i. But then the equations 

Pfc = -£'i£(qfc,qfc+i), Pfc+i = i:'2'C(qfc,qfc+i), 

mean that £(qfc,qfc+i) is just the generating function (of type one, see [GPSOO] for the terminology) of 
the canonical transformation F^. We can then define the right and left discrete Hamiltonian functions 
as generating functions of this canonical transformation of type two and three respectively. Namely, right 
discrete Hamiltonian is 

^"^(qfc,Pfc+i) := Pfc-+iqfc+i - /:(qfc,qfc+i), 
and then the map Fc : T* Q ~> T* Q is given (implicitly) by the right discrete Hamiltonian equations 

qfc+i = L'2'H+(qfc,Pfc+i), Pfc = Di'H+(qfc, Pfc+i). 

Similarly, left discrete Hamiltonian 

■H"(qfc+i,Pfc) := -pfcqfc - £(qfc,qfc+i) 

gives left discrete Hamiltonian equations 

qfc = -i:>i'H"(qfc+i,Pfc), Pfc+i = -i:'2'H"(qfc+i,Pfc). 

For completeness, we mention that the generating function of type four, 

''^(Pfe,Pfe+i) : = qfcPfc - qfc+iPfc+1 + C{qk,qk+i) = qfcPfc - 'H+(qfe,Pfc+i) = -qfc+iPfc+i - "H" (qfc+i, Pfc), 
qfc = L'l7^(pfc,pfc+l), qfc+i = L'2>C(qfc, qfc+i), 

is nothing but the Lagrangian function on the space of momenta satisfying the discrete Euler-Lagrange 
equations 

-D27^(pfc_l, Pfc) + Din{pk,Pk+i) = 0. 
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2.6. Elementary Divisors. For elementary Schlesinger transformations { " } the multiplier matrix Tl{x) 
has a very special simple form, and in this section we describe some useful properties of such matrices that 
we need for our computations. Consider a rational matrix function 

R(a;) = I H , where G — fg^ is a matrix of rank 1. 

X — z 

Because of their role in representing general matrix functions as a product of factors of this form, such 
matrices are sometimes called elementary divisors [Dzh09]. Assuming that tr(G) = g^f ^ 0, which is the 
case we need, we can consider instead of G a corresponding rank-one projector P = f (g^'f)^^g^^, = P. 

Lemma 2.5. Let 

R x) =1+ 

X — Z gTf 

Then we have the following 
(i) basic properties, 

detR(x) = ^, ^[x)-' =1+^/-^ R(z)f ^ f , gtR(^)=('^ 

x — z a; — (, gi f \x ~ z J \x — z 

(ii) the Vanishing Rule, 

(w^R(x)v)f = 0, g^^gt (w^R(a;)v) = 0; 

(Hi) and the Exchange Rule, 

v5v(w^R(x)v) - 9„t(w1'R(a;)v)w^ = Sgt (wtR(a;)v)gt - f5f (w1'R(a;)v). 

Proof. Part (i) is a simple direct computation, see also [Dzh09] . 
Next, notice that for a rank-one projector P = f (g^'f)~-'^g^^, 

f9f (w^Pv) = Pvw^(I - P) and 9gt (w^Pv) g^ = (I - P)vw^P, 

and taking the trace proves the Vanishing Rule. Finally, 

v9v(w^R(x)v) - a„t(w^R(x)v)w^ = vw'fRfa;) - R(a:)vw^ = | ^ ) (vw^P - Pvw^ 

\x — z J ^ 

and 

5gt(wtR(x)v)gt - f9f(wtR(x)v) = f (9gt (wtPv) gt - idf (wtPv)) 



X ~ z 
x~C 



fvwI'P - Pvw^'l 



^x — z ^ 

which proves the Exchange Rule. □ 

3. Schlesinger Transformations 

In this section we derive equations of an elementary Schlesinger transformation { ^ } in terms of the 
coordinates on the decomposition space BxC. We now take R(x) to be of the form considered in Lemma 2.5. 
Then we have the following result. 

Lemma 3.1. Consider transformation (2.2) with 

(3.1) B.{x)=I+^%, z^C. 

x — z g't 

Then 

(i) the poles {ui\ of A{x) coincide with the poles {u,} of A(x) if and only if for some choice of indices 
a, l3 ^ a, fi, and v, we have 

(3.2) z^u^X = U0, f ^ b,,,, ^ b;3,„ gt^<t^c^t^ ^^^ = ^?^:-l, 0-^=0^ + 1, 
where ^ means proportional. Thus, such R(a;) defines the elementary transformation { ^ }. 



(ii) For T{.(x) satisfying (3.2), the residue matrices of A{x) and A(x) are connected by the following equa- 
tions: 

(3.3) AiR(Mi) = R(ui)Ai or equivalently Ai'R.{ui)^^ = 'R.{ui)^^ Ai, i^a,P, 

(3.4) A/sKiu/s) = R.{up)Ai3 and Ac<R(uq)^"^ = R(uq)^^Aq, 
Proof. From the Schlesinger transformation equations 

(3.5) A(a;)R(a;) = R(a;)A(:E) + — ^ or R{x)-^A{x) = A{x)-R{x)-^ 

ax ax 

and 

we see that 

d^^, X 1 / 1 1 \ X / X w X A tr(A,) - tr(Ai) 

tr— R(a;)"^ = = tr A(a;) - tr A(a:) = > 

dx \x — C X — z ) x — Ui 

we see that z, C S Let 2; = Uq and C = up. Then tr(AQ) — tr(AQ,) = —1 and tr(A^) — tr(A^) = 1. 

Since both R(x) is regular and invertible at Ui for i ^ a^P, taking residues of (3.5) at Ui gives (3.3). 
Since R(a;) is regular (but degenerate) at up, we can take a residue at up of the first equation in (3.5) to get 
the second equation in (3.4). Similarly, taking the residue of the second equation in (3.5) at u^ we get the 
second equation in (3.4). 

Further, at Mq, the first equation in (3.5) becomes 

Ua - Up fg^ 




{X - UaY g^f 

^ Up - Ug fgt \ / Ag ^ A^ 

X ~ Ua e^f / \ X — Ua , X — Ui 

The (a; — WQ)~^-terms give 

fgt =fgtA„-A„fgt, 

and so we see that g^ has to be an eigenvector of A^ and f has to be an eigenvector of A^ . Choose an index 
fj, and let g^ ^ c^^^ and f ^ b^^^, then 0^ — 9^ + 1. Taking the residue at u^ gives 

- - A^fg^ _ ^ ^ - fg^A, 

^ Ua-U, gTf ^ Ua-Ui gTf 

but this equation follows from — Aqo = X^ILi ~ SiLi ^ ^Aoo and 

A _A _ '^a-up A.fgt - fgt A, 

which is an immediate consequence of (3.3). Repeating this argument for up and the second equation in 

(3.5) implies that g^' ~ cj^^, f ^ b^^^, and 9p = 9p — 1, for some choice of the index v, which completes the 
proof. □ 

We can extend elementary Schlesinger transformations to the decomposition space as follows. 

Corollary 3.2. Let A{x) and A{x) be connected by an elementary Schlesinger transformation { " } given 
by R(2;) of the form (3.1) satisfying (3.2). Then the corresponding points (B,C^^) and (B,C^) in the 
decomposition space are related as follows: 

(i) If i ^ o:,f3, R(Mi) is invertible and 

(3.6) C| - C^R(mj) { equivalently, Cl ^ ClUiu^y^) and - R(wi)Bi (Bj - R(wj)"^B.,). 
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(ii) For the special indices a and (3 we have 
(3.7) Cl^ClHiup), B^~R(w^)B^, and ^ Cj; RK)-\ B„ ^ RK)-1B„. 

Proof. Multiplying first equation in (3.3) on the right by B^ gives AiR(wi)Bi = R(wi)AiBi = R(ui)Bi0i, 
and so the matrix R(Mj)B,j is the matrix of eigenvectors of A^, R(ui)Bi ^ Bj. The remaining equations are 
proved similarly. □ 

Remark 3.3. Recall that the matrices R(u^) and R(mq)^^ are degenerate, R(u^)b^,i, — 'R.[ua)~^^a,n = 
and c^^R(u^) — c'^'R,{ua)^^ = 0, see (i), so the proportionality coefficients in (3.7) at b^^^, hp^y, cj^^, and 
c^^^ are zero. 

4. The Generating Function 

The main objective of this section is to obtain an explicit formula for the discrete Hamiltonian function 
'H^(B, C^) that generates Schlesinger transformations dynamic on the decomposition space. In this case the 
equations (2.5) take the form (recall Remark 2.3) 



b.. = ||^(B,Ct) and cf.|^(B,Ct), 



where we use the identifications r*(C''') ~ C' x (C^*)^- 
First we introduce some notation. Let 

I = <t<n,l<j< r,} and I' = I\{{a, fi), {f3, j^)}, 

and for a multi-index I — (i, j) £ I' , we put 



(4.1) hj = b,j, 4 = , Oj = ei, Rj = R(u,) when i 7^ a and R„ = I ^ 

Theorem 4.1. Let 



c'^h c^^b« 



X - Ua C^Jh 



(4.2) R(a;)=I- 

(4.3) 5(B, C^; @) = ^ 0i log(4R/b7), and 

leT 

(4.4) n+{B, &; 0) = (0;^ + log(c;^b^,,) + 9^ log(c;^b„,^) + {6^^ - 1) log(cj:tb0_,) + S{B, &; &). 

Then, ifIL{x) is the matrix defining the elementary transformation {'^t} ffom A.{x) to A{x), (B,C^^) S 
[B X C)(0 A^o); (B, C^) E {B X C)(-@ fl'^G tf^^ corresponding points in the decomposition space B x C, and 
&i = CjBi, we have 

(4.5) b,,, ^^(B,Ct.0) cf = ^(B,Ct;0). 

We prove this theorem using a sequence of Lemmas. The generic situation is very simple. 

Lemma 4.2. Let R(x) be any matrix of the form (3.1) with Ui 7^ and let AiR(wj) = R(ui)Ai, with 
A, = B,Cl A, = B,Cl and C^B, ^ cJb, = 0,. Then 

J - d-,, (^91 log (cf R(u,)b.,,) ) and cf = 9b,„ {of log (cf R(?/,)b,., 
Proof. Indeed, in view of Corollary 3.2, C] ~ C|^R(ui), and the normalization Cj^B^ = 0^ gives 

cf n .^^ = d^.. log (cf R(«,)b. 
cfR(ui)bij ^ ^ 

The other equation is proved in exactly the same way. □ 
The following Corollary is then immediate. 
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Corollary 4.3. Equations (4-5) in Theorem J^.l hold for all indices {i,j) with i ^ a, (3 and for {13, j) for 
i ^ V . Moreover, for i ^ a, (3 we have 

A, = X^K.ab,., (H+(B,Ct.0)) ^ A,=J2 (5,.tH+(B,Ct;0)) cf . 

For the index {f3,i'), since R(M/3)by_^ = and c^^R(u^) = 0, a different approach is needed to find 
expressions for hp^^, and c^^. We do it in Lemma 4.5. 

Let us now consider the indices a and f3 defining the transformation. In tliis case R(a;) becomes either 
singular or degenerate and some regularization is necessary. We do it using the equations that cut out the 
locus Ac B xC. 

Lemma 4.4. In the setting of Theorem 4-i with Rq, given by (4-1), we have 

(4.6) b„,^ = a^.t((0^:-i)iog(cj:tb^,,)), 

(4.7) c^^t=aj,. (^Miog^c^tb^ ,^))^ 
and for j ^ 

(4.8) b„,, = {ei log (ci^K^h^^^)) 

(4.9) cit=a^^^^(0Jlog(c^tR^b„,,)). 

Proof. From (3.2), b^^.^ ~ ^i3,u- This and the normalization c^'^b^j,^ = — ~ ^ gives (4.6), (4.7) follows 
similarly from c^^^ ^ c^^. This also gives an additional orthogonality conditions 

(4.10) c^J^fs.u = c^^baj = for j ^ n, and also c^^hp^,, ^ 0, c^^ba,^ ^ 0, 
so Rq is well-defined. 

We now need to solve the equations ^ CJ,R(uq)^^ for Cj, and B^, ^ R(ua)~^BQ, for Bq,. The matrix 
R(wq)~^ = I — b^^^(cj|^^b^ ,y)~^cj^^ has rank m — 1 with the null-spaces spanned by c^^ and Sq.^. We can 
invert it on the orthogonal complements of those subspaces using Rq,. Indeed, 

RqRK)-i=I-^^^ and RK)-iR„-I-^;^. 

Co b^ ,y Cp ba.^ 

Thus, for j ^ a, 'Ra'R-{ua)^^i>a.j ~ i>a,j ^ Raba j, which gives (4.8); (4.9) is similar. □ 

Lemma 4.5. In the setting of Theorem 4-1 with ■H^(B, C^; ©) given by (4-4)7 have 

(4.11) hp^,^d-.,H+{B,&-&), c^^ =aKH+(B,Ct;0) 



Proof. Using (4.3) and the residue at infinity condition — A^o — Tll^i ~ E"=i ~ ~Aoo, we get 
Using the Exchange Rule (iii), whose adaptation for R^ is straightforward, this becomes 
which in turn can be rewritten as 



Thus, 



Using the Vanishing Rule (ii) and the normahzation, we get 



Thus, 



A = At 



which completes the proof of the Lemma (and of Theorem 4.1), and also gives the expressions for hp ^ and 

Remark 4.6. Even though the generating function ?^+(B,C^;0) is defined on a dense open subspace of 
(B. C^), in view of the orthogonality conditions (4.10), to get the well-defined map from {B x C)(@,a^) to 
(B X C)(g)^Aoe) need to restrict it to the subspace 

{B X C^r' = {(B, Ct) I c^hp^, = c^^b„,, = for j ^ fi}. 
5. Examples 

In this section we consider, using our discrete Hamiltonian framework, two explicit examples of discrete 
Schlesinger transformations that can be reduced to discrete Painleve equations. We also try to illustrate in 
this section the role played by the geometry of Painleve equations in not only determining the type of the 
equation, but also in studying the relationship between different explicit forms of equations of the same type. 

Recall that, according to the geometric theory of Painleve equations developed in [SakOl], both continuous 
and discrete Painleve equations can be classified by the type of their space of initial conditions, introduced 
by Okamoto [Oka79] in the continuous case. This space can be obtained by successively blowing up the 
projective plane at nine (possibly infinitely-close) points that lie on a unique (possibly singular) cubic 
curve. The type of the resulting surface X is then given by the configurations of irreducible components 
Di of its unique anti-canonical divisor D or, equivalently, by the root sub-lattice R spanned by Di in the 
orthogonal complement of D in the Picard lattice Pic(X). The orthogonal complement R-^ of R in 
D-^ then corresponds to the symmetries of the equation — its Weyl group acts as automorphisms, called 
the Cremona transformations, of Pic(Ar) preserving R C and the translational part of the Weyl group 
corresponds to discrete Painleve equations. 

We also want to emphasise the following point, see also [Sak07]. Given a discrete Painleve equation, 
we can determine its type by constructing its space of initial conditions. However, there is no converse 
process, although there are some particularly nice forms (that we call standard) of equations of each type, 
see [Sak07, GRO03, Mur04]. Further, there can be many non- equivalent equations (i.e., equations that can 
not be transformed into each other by some change of coordinates) of the same type — such equations can 
correspond to different translation directions in the affine Weyl group of symmetries of X. And even when 
the two different equations are equivalent, finding the actual change of coordinates that will transform one 
equation into the other can be a very challenging problem. In each case, understanding the geometry and 
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the blow-up structure of the space of initial conditions can be of great help in determining the relationship 
between two different equations of the same type, as we illustrate by the two explicit examples below. 

For the first example, we consider the Fuchsian equation of matrix size m = 2 with n = 3 finite poles and 
rank-one condition ri = r2 — — 1 for the residue matrices at those poles. It is well-known that continuous 
isomonodromic deformations of this system are described by the classical continuous Painleve equation Pyi, 
the type of the space of initial conditions is and the discrete Painleve V equation d-Py = d-P{D^'^) 

appears as Backlund transformations of Pyi- Thus, it is reasonable to expect that an elementary Schlesinger 
transformation of this system is described by d-Py, and we explicitly show that this is indeed the case. 

The second example is a more interesting one. In [Sak07] one of the authors (H.S.) posed a problem of 
writing down discrete Painleve equations of type A^'^** , A^^"^* and A'^'^* as Schlesinger transformations of 
some Fuchsian equations. One difficulty here is that these equations do not correspond to any continuous 
Painleve differential equations. This problem was solved theoretically by Boalch in [Boa09], where it was 
shown that the moduli spaces of monodromy groups of some Fuchsian equations are analytically isomorphic 
to the spaces of initial conditions for these discrete Painleve equations. In particular, Boalch showed that the 
Fuchsian system with matrix size m = 3 with n — 2 finite poles (and so there are no continuous deformation 
parameters) and rank-two residue matrices, ri ~ r2 ~ 2, has as its symmetry the afhne Weyl group of type 
Eq^\ and so the type of the corresponding surface is either ^2^^ (for the multiplicative or q-dijference case) 
or ^2^'* (for the additive or difference case). Since we consider the additive case, we expect that discrete 
Schlesinger transformations of this system should have type d-P(A2^^*). However, an elementary discrete 
Schlesinger transformation, when written in coordinates, becomes very complicated and in fact it can not be 
reduced to a standard form of d-P{A'^'^*), as written in [Sak07]. Using the geometry of the space of initial 
conditions, we see that this is an example of a situation when the standard form of the equation corresponds 
to a direction that is a combination of two different directions given by different elementary Schlesinger 
transformations, and then we explicitly show how to combine them into a single Schlesinger transformation 
that corresponds to the standard form of d-P{A':^^*), this is a new result. The geometric approach that 
we use requires finding correspondences between the rational surfaces that appear as desingularizations of 
different coordinate forms of two difference equations of the same type. Such correspondence gives us a blow- 
down structure that can then be used to find a change of coordinates that relates the two equations that we 
consider, see the recent preprint [CT12] by A. Carstea and T. Takenawa for a more detailed introduction 
into this technique. 

5.1. Discrete Painleve V equation (type £'4^''). 

5.1.1. Elementary Schlesinger transformations. For this example we consider a Fuchsian system of the spec- 
tral type 11, 11, 11, 11, i.e., we have n = 3 (finite) poles and the matrix size is m = 2. In view of Assump- 
tion 2.2, it is possible to make the rank of the residue matrices at the finite poles equal to one. Further, 
using the Mobius transformations on P^, we can map (uq = 00, ui, U3) to (00, 1, 0). We put U2 ^ t and from 
now on we use (1, t, 0) instead of the indices (1, 2, 3), so O2 — 9t, etc. Then 



A(x) = ^ + ^ + ^, Ao, = -(Ao+Ai+Ai 

X — 1 X — t X 



Kl 
K2 



and we have the following Riemann scheme and the corresponding Fucks relation: 
X = X = 1 X — t X ~ 00 ^ 

0Q Ol Bt Kl \ , 00 + 01+61+ Kl+ K2= 0. 

K2 } 

In the rank-one case there is no need for the second index, and so we put = b^.i, cj = c|^. Then the 
spectral decomposition is simply A^ = b^cj with tr(Ai) — c|bi = 0i. 
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Consider the elementary transformation { J *}, i.e., 0i = 9i — 1, 9t — Ot + 1 and da = 0o. Its multiplier 
matrix and the generating Hamiltonian function are 



•H+(B,C;0) 



1 - t htcl 



= log 



+ Oi log 



• (clbi)(clbt) ' 



J + 9t log (clbt) + 00 log (4R(0)b( 
1) log (c^bt) + 0, log (c|bi) 



6*0 log (^c|bt cjbo - (1 - t)c|bo cjb; 



This function generates the following dynamic equations: 



(5.1) 
(5.2) 
(5.3) 
(5.4) 
(5.5) 
(5.6) 



bi 
bo 



d_^,H+^{0,-l)^, 
' c^bt 

{clht)ho - {1 ~ t){c%)ht 

'(clb0(4bo)-(l-t)(clbo)(cjb0 

b* , . bi _ 



+ _ 



CA = 



= "70 



l-Oi- 



R(0)bo 
4R(0)bo' 
(4bo)bt - (1 



t)(4bt)bo 



c b 



cjbi (cIbt)(4bo) - (1 - t)(cIbo)(cSbt) ' 



C bi 



= Co 



(clb,)4-(l-<)(4bt)cl 



(clbt)(4bo) 



+ 1 - 6ii - 00 



(l-0(4bo)(4bt) °4R(0)bo 



4r(o) 



clh 



t "t 



(4bo)4-(i-t)(4bo)4 
cib, ' '''^(4b,)(4bo)-(i-t)(4bo)(4b, 



Remark 5.1. To illustrate the general situation, we now show how to solve these equations to define a map 
from {B X C)(0.Aoo) to {B x C)(@ ^^-j in this relatively simple example. This map is defined up to trivial 
transformations and becomes unique once we specify some normalization condition. 

Indeed, equation (5.4) implies that c| = Ctc\, where Ct is some proportionality constant. Thus, we can 

rewrite the multiplier matrix as R(a;) = I + -^-zj-jr^- Equation (5.5) can be rewritten as cj ^ CqR(O), and 

so Cq — co4R(0)"\ cf. (3.6). Finally, expression for c[ can either be obtained by a direct computation or, 
simpler, from the equation c|Aoo = c|Aoo. We then substitute these expressions into equations (5.1-5.3) 
and simplify. The final map is given by the following set of equations (with cq, ci, and ct arbitrary constants 
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corresponding to trivial transformations): 





bfC 



c b 



Let us now look at the Schlesinger dynamic on the space of accessory parameters. Normalizing the 
eigenvectors bj and c| , we can write 



A, 



[c. 1] 



w 



[j, +w 1] , tr(A, 



a-iici + bi) = ai{/3i + ji). 



Here w — —bo and so /3o = and 70 = ^o/^o- The reason for introducing w is to simplify the relations between 
parameters Oi, bi, and Ci that follow from the trace conditions and the equation Aqo = — diag{/ci, ^2}- These 
relations, written in terms of a^, /J^, and 7^, are 



(5.7) 
(5.8) 
(5.9) 
(5.10) 



ao + ai + at = 
aiA7i + atl^tlt + w{ki - ^3) = 
aiPi + atPt + K2 = 
6*0 + ai7i + o^tlt + Ki = 0, 



where the last equation also follows from (5.9) and the Fuchs relation. Thus, everything is parameterized 
by oi, Pi, at, and /3( subject to one additional constraint ai(3i + at(3t + K2 = 0. One additional degree of 
freedom is related to the global gauge action by constant non-degenerate diagonal matrices that allows us 
to normalize one of the a^s, but it is convenient to keep this freedom for computations. Finally, note that 
the normalization parameter ai should be considered as a part of either b^ or cj, where the actual choice 
depends on the context. 

Same coordinates were used in [SaklO] to describe the isomonodromic dynamic in the continuous case 
and the reduction to Pvi ■ With ai attached to c\ , the standard symplectic form on the decomposition space 
reduces as follows, see Proposition 4 in [SaklO] for a more general statement: 

3 

w = V"tr [dC] A dBi) = dao A d{-w) + dai A d{(5i - lu) + dat A d{l3t - w) 



1=1 

d{ao + ai + at) A d{—w) + dai A c?/?i + dat A d/St = dai A c?/3i + dat A d/St 

'dat dai 
at ai 



dai ,/ r, ^ dat 
— A d(ai/3i) + — A diatPt) 
ai at 



A d{atl3t 



Ad{ail3t), 



where we used (5.9) to eliminate aiPi. It is convenient to introduce new symplectic coordinates p ~ (ai/3t)/t 
and q — —{tat)/ai. This gives, using (5.7-5.10), 

ai(ii=pq-K2, ail3t=tp, aiji ^ 9i + K2 - pq, 

atfit^-pq, atl3i ^ '^{k2 - pq), atjt^6t+pq- 



15 



For the continuous deformations, the only deformation parameter is t and the corresponding Hamiltonian, 
given by equation (2.4), under this reduction becomes 

tr(AtAo) , ti{AtAi) _ {Ot + pq){eo + p{q - t)) {Oi + K2 - p{q - t)){tet + K2q - pq{q - t)) 



Ht = 
Since 



t 



t 



tit -I) 



A d{tp) ^ dpAdq-d (Jj^ A dt, 



uj-dnt Adt = dpAdq-d['H 



pq 
t 



Adt = dp Adq- dHvi A dt, 



where Hvi = T~L 



-pq/t is the Hamiltonian function for Pyi- 
Unfortunately, such a reduction procedure does not work in the discrete case and we need to work directly 
with the evolution equations. We give a brief outline of such computations and show the reduction to d-Py . 

It is convenient to attach the normalization parameters di to b.^ and to c| (i.e., in what follows, b.^ 
denotes the normalized vector, but in the formulas describing the dynamics we use Oibi, and so on). It turns 
out that it is possible to separate the "normalization" dynamics of a.i and the "directional dynamics" of /3i 
and 7i, the latter one is given by the Painleve equations. 
Equation (5.1) becomes 



(5.11) 



aibi — fli 



1 




1 


A - ^. 




Pt - w_ 



clb, 



and so 



w = Pi- (3f 



Note that c}bt = 71 + w + /3t 
Equation (5.4) becomes 



■ w 



7i + /?! = c|bi, and so 6*1 — 1 = aicjbt — aic|bi = 9i, as it should be. 



(5.12) 



ai 7i 



and so w 



71-7* 



Using this equation in the computation of the trace, c|bi — ji + Pi = c\hi, we again see that the trace 
condition is satisfied, aic|bi = di. Important consequences of (5.11) and (5.12) are 

(5.13) clbt=7i+ft=7t+/?i=clbi, 

(cjbt)(w -w) = {^t+w + Pt- w){w -w) = 7i(/3i 

We show shortly that the first equation generates the "normalization" dynamics 
Equation (5.2) becomes 

1 1 ^0 



i) + A (71 - It) = /3i7i - Pat- 



ao 



-w 



(l-i)(cIbo)(cJbt)-(cIb,)(4bo) 
The traces in this formula are, using (5.11) and (5.12), 



(I-O(clbo) 



cJbo 



cjbt 
cjbt 
Cobo 



It - 

70 ■ 

71 ■ 
70 ■ 



71, 



■Pt 



w 



70 + Pi 



70 



1 

Pt-w 



- aoPi 



(clht) 



ao 



c% 



= 7t + /3i 
00 + do{w - 



= ctbi, 

w) 



ao 



The equality of the normalized vectors then gives 

(1 - t)(7i)(-w)) - (71 + Pt){~-w) - (1 - 0(7i)(A -w)- (71 + Pt){-w), 
(1 - t)(7i)(-^i) = -(71 + Pt){w -w)^ Pat - 7iPt, 
which gives the following important relation 

(5.14) Pat = tpai ■ 

As before, the normalization condition reduces to — Bo- 
rn 



It is possible to show that the remaining evolution equations do not give any additional information. 
Let us now write everything in the symplectic coordinates p and q. Equation (5.13) gives the equation 
for the normalization dynamic, 

(k _ Qt(7t + A) _ Ot+Pq+ 1(^2 -pq) 

ai ai(7i+/30 di+K2-p{q-t) 
Multiplying (5.14) by aidt, we get 
(5.15) (6*1 + K2~ pq)q{fi2 - pq) = tp{6t + pq). 

Finally, combining (5.11) and (5.8) gives 

ai/3i(7i + A) + Ptiat-it - ai/3i + ki - K2) = ai^t(7t + ,^1) + h{aili - atPt + ni ~ K2). 
Dividing by 71 + /?< = 7t + /3i and doing some rc-arranging, we get 

(5-16) pq+pq~K2^ ^ e,+., + ^ e,+., - 

We now claim that after a certain change of variables equations (5.15, 5.16) transform into the standard 
form for d-Py . Although this can be seen directly without too much difficulty, we prefer to use the geometry 
of the rational surface corresponding to the equation to explain how to find this coordinate transformation, 
illustrating some of the ideas that we need for the next, more complicated, example. 

5.1.2. Geometry of d-P{D^^^). Let us briefly review the geometric description of the difference Painleve V 
equation from the point of view of the theory of rational surfaces [SakOl]. We think about it as a birational 
map : pi X pi pi X pi with parameters oq, . . . , 04, s or, equivalently, as an automorphism of a field 
C(ao,...,a4;s;/,5), 

ai a2 03 ^ \ fai 02 03 _ j _\ ai=ai+5, a2 = 02 - S, 03 = 03, 

; s; /, g 1— > _ _ ; s; /, g , _ _ j. 

04 flo J \ 04 flo J a4 — 04, Go = Oq + 0, s = s, 

where S — ai + 2a2 + 03 + 04 + 00 and where / and g are given by the equation that is usually called the 
discrete Painleve-V equation 

ai , ao 



(5.17) 



99 



g + 1 sg + 1 
(/ + a2)(/ + a2 +04) 



sf{f - 03) 

we refer to it as the standard form of d-Py- We lift this map to an isomorphism (p between rational surfaces 
Xa — X^ao,...ai) and Xa obtained from P^ x P^ via a sequence of blowups that resolve the indeterminacies 
of the map Lp, 

: *" 

I f 
I 

Y ¥ 



X r" 3- r" X 



Decomposing 1^9 as a sequence of two maps, (/, g) (/, g) (/, g) it is easy to see that the indeterminate 
points oiipi are (00,— l)and ((X),— l/s)and the indeterminate points of(/52are (—02,0), (—02 — 04,0), (0,00), 
and (03,00). Applying (p^ we obtain that the indeterminate points of (/? are pi(— 02 — 04,0), p2(— 02,0), 
7)3(0,00), p4(a3,oo), 7)5(00,-1), and P6(oo,— 1/s). In exactly the same way we find that the indeterminate 
points of ip^^ are j)i(— 02 — 04, 0), . . . ,pq{oo, — 1/s). 

As usual, we denote a local coordinate description of a blowup at a point {xo, yo) by (x, y) -s— (u, v) U {U, V) 
where x — Xq ^ u = UV and y ^ yQ — uv — V , then the exceptional divisor is given by u = = 0. For 
example, extending the map ip at the blowup at pi (— 02 — 04, 0), we use ui =7 + 02 + 04 and uivi = g. Then 
/ = f{ui,vi) and g = g{ui^vi) become well-defined and we see that the exceptional divisor Ei is mapped 
to the line / = —02. On the level of divisor classes on Pic(X) we get the map Ei H> Hf — E2, where Hf is 
the class of the total transform of the line / = const. 
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Proceeding along the same lines and notieing that we need to blow-up at two additional points ^7(0,01) 
on E5 and ps{0,ao/s) on (here we give the coordinates of the points in the {ui,Vi) charts), we get the 
following blowup diagram describing the construction of the surface X^, corresponding to our equation: 




The Picard lattice of is generated by the total transforms Hf and Hg of the coordinate lines and the 
classes of the exceptional divisors Ei, 

8 

Pic(X) = ZHf e Zi?g ® ZEi. 
The anti-canonical divisor —Kx uniquely decomposes as a positive linear combination of — 2-curves Di, 

8 

-Kx = 2Hf + 2Hg-Y,Ei = Do + Di+ 2£»2 +D3 + D4, 

where 

-Do = Hg — El — E2, Di = Hg — S3 — S4, D2 = Hf — E5 — Eq, Dz = E5 — E^, D4 = Eq — Eg. 

The configuration of these —2-curves (on which the blow-up points are located) is described by the affine 
Dynkin diagram of type 




which is why we say that the difference Painleve V equation has the type d-P{D\^ ). The corresponding 
divisors Di generate a sub-lattice R = Span^lDo, . . . , D4} C (— iCx)"'" — E^^^ in Pic(X). Symmetries of 
X are given by the affine Weyl group of Cremona transformations generated by the complementary root 
sub-lattice in {—Kx)'^ acting as automorphisms of Pic(X). In this case, R^ = Span^jao, . . . , Q!4}, where 

ao = El — E2 
cxi = E^ — E4 

a2=Hf-Ei- E3 

Q!3 Hg — Eq — Ej 

0:4 = Hg — Eq — Eg 
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which is again of the type . FinaUy, we can compute the action of on Pic{X) to be 





^ 5Hf - 


V2H,- 


- 2Ei 


-2E2- 


2E3- 


2i?4 — E^ — Eq — Ej ~ Eg, 


El ^ Hf 


-E2 




^ 2Hf - 


hHg- 


El - 


E2 ~ E^ 


— E4, 




E2 1^ Hf 


~Ei 


E5 


^ 2Hf - 


hH,- 


El - 


E2 — E^ 


-Ei- 


- Es, 


Ez ^ Hf 




Eq 


^ 2Hf - 




El - 


E2 — E3 


~Ei- 


-Er, 


E4 Hf 


-E3 


Ej 


^2Hf- 




El - 


E2 — E^ 


^E^- 


- Eq, 






Eg 


^2Hf- 




El - 


E2 — E3 


-Ei- 


-E5, 







and so the induced action ^p.^, on the sub-lattice R is given by translation: 

ai I— > ai as as — 6 6 = —Kx- 



5.1.3. Reduction to the standard form. We now demonstrate how to reduce equations (5.15, 5.16) to equations 
(5.17). First, we resolve the indeterminacies of V' '■ (PiQ) ~^ (PiQ) to get an isomorphism tp : Xa — > Xa- 
The surface X — X^, with a denoting the collection of parameters of the map, o — {Sg, ^1, ^t, '*2}j is 
constructed by the following blow-up diagram: 



p = 9 p = 00 



P3(0, 00) 



P6(0O,l) 

pi(0, 1/(01 + ^2)) 

P5(0O,t) 4 



pi(oo,0) 



I pai^i, 1 - Ki + K2) 
■* — 



ip7(0,-6'o) 



g = 



p1{Q--h) 




Here the coordinates of the blow-up points are given w.r.t. the following local coordinate systems: 



Pi I ^ - 0,9 = 



P3 ( P = 0, - = 



P5 - = 0,9 = t 



P6 - =0,9 = 1 



- = 0,_p(7 
1 



P8 



Jt \ and — = 0, 9 = 



P4 ( P = 0, — = — ) and P4 (p = 0, — 
pq K2J \ pq 



O,piq-t) = -eo)j , 
0,p{q - 1) = 1 - Ki + K2) 



O1 + K2J ' 



Note that the resulting configuration of exceptional divisors looks quite different from the previous case 
and moreover, rank(Pic(X)) = 12. The reason for this is that the surface X is not relatively minimal for ip, 
see also [CT12] — there are two pairs of — 1-curves that are exchanged by tp: F2 O F^ and F2 O Hq — F^, and 
so blowing down any of these pairs will still result in an isomorphism between the corresponding surfaces. 
We blow down F2 and Hq — F3, which eliminates —4 and — 3-curves and gives us the diagram on the left. 
This diagram now looks very similar to the configuration we had for the standard form of d-Py (on the 
right), 
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and in particular we immediately see that X is again of type D\ '. We now need to find the blowing-down 
structure for X to match the blowing-down structure for X. Thus, in Pic{X) we put £g = F2, £10 — Hg — 
and we are looking for rational 0-curves Hf, Hg, and — 1-curves £i, i = 1, . . . ,8 such that 

UfHg^l, £f^-l, 'H) = 'Hl = 'Hf£i='Hf£r^£i*£j^Q, 1 < « ^ j < 10. 

Using the above diagrams, the genus formula 2g{C) — 2 = + C • Kj^, and the distinguished — 2-curve D2, 
we see that it makes sense to take Hf = Hp + Hg — Fi — and £i — Fi, i = 5, ... 8. Then the possible 
choices for Hg are Hg, Hp + Hg ~ F2 — F3, Hp + Hg — F4 — F3, and Hp + Hg — F!^ — F^. Attempting to 
represent Fi — F2 = Hg — £1 — £2 eliminates the first two possibilities, and the remaining two are equivalent. 
To summarize, we have the following blowing-down structure: 

^/ = Hp + Hg — Fi — F3, £3 = Hp — F3, £i — Fi, i — 2,4, ... ,8, 

Hg = Hp + iJq — F3 — F4, £q = F2, 

£1 = Hp + Hg — Fi — F3 ~ F4, f 10 = Hg — F3. 

To find the expressions for coordinates / and g we proceed as follows. The linear system \Hf \ consists 
of curves of bi-degree (1,1) passing through pi an pa, and so \Hf \ = {cipq + C2 ~ 0}. Thus, [pq : 1] is 
a projective coordinate on \Hf\, and so we can take / = pq. Similarly, [pq — (^1 + K2) '■ p] is a projective 

coordinate on \Hg\. We take g = —j (^q — ^^^^^ , where the additional factor — 1/i is needed to ensure that 

the (/, (7)-coordinates of are (00, —1), as in the standard case. 

We can now match the parameters of our elementary Schlesinger transformations with the parameters s 
and fli in the standard form (5.17) of d-Py. For example, £1 blows down to (/, 5) = (^1+^2,0) = (—02—04,0). 
Doing this for other divisors £i results in the following 

ao — 61 - 1 + ki, ai — -~9t - fci, 02 = Of, as = K2, 04 = -9i - 9t — K2, ^ = -1- 

Thus, not only we can now immediately see that equations (5.15, 5.16) transform into the standard form 
(5.17) of d-Py, but also that under this identification of parameters, 

9i = —0,2 — as — 04 = —02 + (5 — 03 — 04 = 01 — 1, 9t = —So — oi — 02 — 03 — 04 — 1 = 0t -f 1, 

and the other parameters remain unchanged. Thus, we can identify the action of d-Py on the Riemann 
Scheme of the Fuchsian system under this correspondence. We use the same observation in the next example. 

5.2. Discrete Painleve equation A-P{A''2^*). 

5.2.1. Schlesinger Transformations. We now consider Fuchsian systems of the spectral type 111, 111, 111, 
which corresponds to n = 2 poles, matrix size m = 3, and, rank(Ai) — 2. Thus, 

r„it" 

A, = B,C^ = [b,,i b,,2] 2t • 
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Using a Mobius transformation preserving u = oo, we can map the poles to to Ui = and U2 = 1, then 
the Riemann scheme and the Fuchs relation are 



el 
el 



K3 



^\ + el + 61 + el 



2^ ^3 



This example does not have any continuous deformation parameters but it admits non-trivial Schlesinger 
transformation. As before, consider an elementary Schlesinger transformation { \ \ }, e\ =6\ — l, 62 = ^2 + !> 
and 9l = Oj otherwise. 

The multiplier matrix R(a;) and other related matrices in this case are 



R(a;) = I 



Ri = 1 



1 b2,iC2^ 

^C2^b2,l^ 
bl,lC2^ _ 
C2^bi,l 



b2,ic}^ 

.It 



C,'h2 



R(a;)-^ =1 + 



b2,lC2''' 



a;-lc^^b2,i 
R2 = R(l) = I - 



b2,lC2^ 
C2^b2,i 



R(0)-\ 



and we immediately see that 



If 



R1R2 — I 



b2ac^ 

Ci^b2,l ' 



R2R1 — I — 



bi.ic}'^ 
Ci^bi^i 



The generating Hamiltonian function can be written as 



n+{B, &; 0) = - 6[ -0'l + i) iog(c^tb2,i) + {el -ef-i) iog(c}^b2,i) + el logPi + ei iog2?2, 



where 



= (C}^b2,l)(c?tbi,2)(c^^bi,i) - (C}^b2,l)(cf bi,i)(cl^bi,2) " (c}^bi,2 ) (cf b2,l) (c^^bi.i ) , 
T>2 = (C2'''b2,l)(C2'''b2,2) - (C2^b2,2) (C2^b2,i ) . 



We restrict this function to the subspace {B x given by the equations c^^b2,i = C2^bi,2 = 0. It then 

generates the map from {B x C)(0^Aoo) *o ^ 0(6, Aoo)' where Aqo = A^o = diag{K;i, k,2, ^3}. The map is 
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implicitly given by the equations 



iol-i) 




, a2 (C2^b2,2)b2,l - (C2^b2,l)b2,2 

+ V, ' 

2 (C2'''b2,l)b2,2 - (C2^b2,2)b2,l _ 2 R'2b2,2 



V2 'cVR^h^ 



202,2 



^1 




i C/bi,2 

c}'''b2,i bi,2 



,2f _ ^2 (^2 ^2,1)02^ - (C2^b2,l)c2^ _ 

1^2 C2'R2b2,2 



cf. (3.6-3.7). These equations can be solved to give the map from [B x C)(0.a^) to {B x C)(@ ^.^j explicitly, 



(5.18) C}t = c\ [ C^t + (^1 _ ,1 _ l)_£^ + 4!^cr + ^l^(b2,lCr+b2,2Cr)b,2 ^2tj, 



(5.19) 


_ „2 2t 
^1 — '^I'-l 


(5.20) 


'^2 — '-^2'^1 1 


(5.21) 


^2^ = C2 I C. 



b2,lCi 



If 



c?cf R2, 



2t 



c}^b2,i^ 

C2^(bi.iC^'^ + bi.2cf )b2.1 c}^ 

{ei-e\-i) c}tb2,i^ 
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(5.22) 
(5.23) 
(5.24) 
(5.25) 



bi.i = 
bi,2 
b2,i = 

bo. 2 = 



3-b2,i, 



cJ^(b2,lC2'^ + b2,2cf )bi,2, 

"1,2 — . i-\ ^ — '-*2:1 



bi. 




\ . .^ b2.1 



Ci b2,l 

Ci^b2,l cj^b2,i 



bi, 



b2,2 — ^R'2b2.2- 



C2^(bi,ici^ + bi^2cf )b2,l 



2t^ 



(?1-I)(c}^b2,l) 



R2b2 



We claim that the resulting dynamical system has the type A2 , as expected. To show this we first need 
to introduce some coordinates on the phase space. Unfortunately, using symplectic coordinates similar to the 
previous example results in very complicated equations, so we use the following parameterization instead. 

Assuming the generic situation, we can use the similarity and trivial transformations, see Remark 2.4, 
to first map the vectors bi^i, bi_2, and b2,i to the standard basis vectors using the global similarity trans- 
formations, and then use the trivial transformations to adjust the scale on the coordinate axes so that all 



components of b2_2 are equal to 1. Then, using the condition Cj^B^ 
zation: 



©i, we get the following parameteri- 



Bi 



1 
1 




6\ 
0? 



B2 



1 

1 

1 1 



-y y + 



where the discrepancy in the notation of the elements of C\ and C\ is due to the fact that we can express 
a and /? in terms of x, y, and the parameters 9l and Hi by comparing the coefficients of the characteristic 
polynomial of Aoo = — BiC| — B2C2. The resulting expressions, although easy to obtain, are quite large 
and we omit them. Using equations (5.18-5.25) we get the following dynamic in the coordinates {x,y): 



a 



V = 



{a - p)y - m 



1) 



1 



(/3 - a)x + {ef + l + a){e^ + a) 
{a-l3){y + ej)-a{e'2 + l) 



where we still need to substitute a = a{x, y) and /3 = l3{x, y). 

With the help of Mathematica we can find the indeterminate points of the map -0 '■ {x,y) ~> {x,y) to 
get the following blowup diagram (note that in this case it is slightly more convenient to compactify to 
p2 rather than x P^), 



(5.26) 



B:z = 0. ,E :x = 




Q = 2E 




E — E2 — Ei 



E — Eg — Ej 
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where the projective coordinates of the blowup points are 

(1 : 1 : 0) , P. {01 : : 1) , ^^^^''^ + + ■ + '^^^^^^ + + 

pe (0 : -ei : 1) , p7 (0 : 1 : 0) , p4 



9? - 01 



n 01 



n " d\ 



n - o\ 



n-o\ 



'y-x Z 



n - d\ 



P8 (1 : : 0) , P9 0, ) , P5 I : : 1 



Pg is given in the coordinates (u, v) = ' y-x ) ' ^^i*^ ^^"^ quadric Q passing through the points pi, . . .pe 

is given in the homogeneous coordinates [X : Y : Z] hy 

(5.27) {el - el){x -Y- elz){x -y - elz) + {e\ - Bl)z{Bl{e\z - x) + b\y) = o. 

The anti-canonical divisor —Kx uniquely decomposes as a positive linear combination of — 2-curves Di, 

9 

-Kx =?.E-Y,E^ = Do + D^+ D2, 



where 



Dq — 2E — Ei, Di — E — El — Ej — Eg, D2 — Ei — Eg. 



Thus, we see that the configuration of components Di is indeed described by the Dynkin diagram of type A^2^ . 
To this diagram correspond two different types of surfaces, the generic one corresponds to the multiplicative 
system of type ^2^'*, and the degenerate configuration where all three components Di intersect at one point 
corresponds to the additive system denoted by ^2^^*, this is our case: 
^0 




Dynkin diagram A2 



(1) 





surface 



^2^ ■'-surface. 



Thus, we have showed that an elementary Schlesinger transformation of this Fuchsian system is given by the 
discrete Painleve equation of type d-P{A^2^*). However, equation (5.2.1) is very complicated, and we again 
consider whether it is possible to transform this equation into a standard form. 

5.2.2. Geometry of d-P{A^2^*). Following Murata [Mur04] and Sakai [Sak07], we take the standard form of 



the d-P(A2 ) to be the following. We consider d-P{A2 ) to be a birational map (p : 
with parameters 61 , . . . , fog 
(5.28) 



61 62 
^5 be 



b7 bs'^'^nVh h 



bs = bs, 

64 — &4 



bi = bi, 

bj Bg'^'^j' 62 = 62 
5 = bi + ■ ■ ■ + bg, and / and g are given by the equation 

(g + bi)ig + b2)ig + b3){g + bi) 



^5 = ^5 + S, 
^6 = ^6 + (5, 



b-j — b-j — 5, 
bs^bs- (5, 



(5.29) 



(/ + .9)(/ + .9) 
{f + 9){f + 9) 



{9-b5){g-be) 

(/~^l)(/-&2)(/-b3)(/-&4) 

{.f + b7-S){f + bs-S) 



Note that equation (5.29) is the asymmetric version of equation (1.3) in [GRO03]. We get the following 
blowup diagram describing the construction of the surface Xb corresponding to our equation: 



(5.30) 



P7(-&7,00) 



Ps{-bs,oo) 
• — ^ •■ 




E7 — Ea 



As in Section 5.1.2, Pic(Xb) = ZHf ® ZHg ® 0^^^^ ZEi, The anti-canonical divisor -Kx = 2Hf + 2Hg 
Si=i decomposes as —Kx = Dq + Di + D2, where 



Do — Hf + Hg - El - E2 — E3 - i?4 
The symmetry sub-lattice — Span^jao, . . . ,ae}, where 



Di ~ H f — E^ ~ Eq 



"0 


-i?3 


— i?4, 


ai 


= E2 


-E3, 




a2 


^Ei 


-E2, 


as 




-El 


-Er 


ai 


= Ej 


- Es, 




= Hg 


-El 


-E5 


ag 


= E, 


-Ee 











D2^Hg 

O «6 
6 "5 



E.. 



ao 



-o 



o o 



Qf2 



as 



-O E, 

a4 



(1) 



which is of the type Eq'^K We compute the action of ip^ on Pic(X) to be 



(5.31) 


Hf 6Hj - 


^SHg- 


' 2Ei 


-2E2- 


2E3- 


- 2E4 


-E5 


-Ee 


(5.32) 


Hg ^ 3Hf - 


^Hg- 


Ei- 


E2 — E3 


-E4 


-Er 


— Es, 




(5.33) 


El 1-^ 2Hf - 


^Hg- 


E2- 


E3 — E4 


-E^ 


— Es, 






(5.34) 


E2 ^ 2Hf - 


^Hg- 


Ei- 


E-i — E4 


-Ey 


— Es, 






(5.35) 


Ez ^ 2Hf - 


^Hg- 


Ei- 


E2 — E4 


-E, 


— Es, 






(5.36) 


E4 2Hf - 


^Hg- 


Ei- 


E2 — E3 


-E^ 


— Es, 






(5.37) 


£'5 i-> 3Hf - 


^Hg- 


Ei- 


E2 — E3 


-E4 


-Ee 


-Er 


-Es 


(5.38) 


Ee ^ 3Hf - 


^Hg- 


Ei- 


E2 — E3 


-E4 


-E5 


-Er 


— Es 


(5.39) 


Eri-^ Hf- 


Eg, 














(5.40) 


Es^Hf- 


E7, 















3Er — 3Es, 



and so the induced action (p^ on the sub-lattice is given by the following translation: 

{ao,ai,a2,a3,a4,a5,ae) (ao, ai, 02, as, "4, "5, "e) + (0,0,0, 1,0,-1,0)6, where 6 



-Kx. 



5.2.3. Reduction to the standard form. To match the surface X described by the blow-up diagram (5. 20) 
with the surface Xb described by diagram (5.30), we look for the blow-down structure in Pic(X), i.e., we 
look for rational classes T-lf, Hg, £1, . . . ,£s in Pic(X) such that 



Hf • T~Lg — 1, Sf 



Hi 
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Hg • £i — £i • £j — 0, 1 < i, j < 



and the resulting configuration matches diagram (5.30). Comparing the — 2-curvcs on both diagrams, 

Dq = 2E — El — E2 — E3 — E4 — — Eq = TLf + T-Lg — £1 — £2 — £3 — £4, 
Di = E — El — Ej — Es = Hf — £5 — £&, 
D2 = El — Eg = Tig - £7 — £s, 

we get the following identification (it is easy to see that it satisfies the genus and the rationality conditions) : 

Hf = E — El, £1 = E3, £3 — i?5, £5 = Et, £7 ~ E — El — E2, 
Hg = E — E2, £2 — E4, £4 = Eq, £q = Es, £s, = -Eg. 

To complete the correspondence it remains to define the base coordinates / and g of the linear systems 
\'Hf \ and \'Hg\ that will map the exceptional fibers of the divisors £i to the points tt^ such that tts and TTg are 
on the line / = 00, Try and TTg are on g — 00, and tti, . . . 7r4 are on the line f + g = 0. Since the pencil ["H/j 
consists of all curves on passing through pi(l : 1 : 0), 

= \E-Ei\ = {aX + bY + cZ = 0|a + 6 = 0} = {a{X - Y + PZ) + cZ = 0}, 

we can define the projective base coordinate as /i = [X — Y + PZ : Z], with /3 being some parameter. 
Similarly, 

\ng\ = \E-E2\ = {aX + 6y + cZ = I 6*^0 + c = 0} {a{X - 0lZ + aY) + bY = 0}, 
and gi = [X -e^Z + aY -.Y]. Then 

/i(7r5) = fi{p7) = [-1 : 0], 31(^7) = giiPi) = [1 + a : 1], 
fiM = flips) = [1 : 0], giiiTs) = giipi) = [1 + « : 1]. 

Thus, we take a = —I (and looking at the other points we notice that it is also convenient to choose 
/3 = ^(^2) J so we provisionally put 

X -Y -ejz Y 

h = ^ , 92 



z ' x-Y-e\z' 

In these coordinates equation (5.27) of the quadric Q becomes 

z\f2 + el - e\) {{e\ - el)f2 + {el - el){(e\ - el)g2 - eD) - o, 

and so the points tti, ... 714 lie on the line {6\ — el)f2 + {e^ — ^2)((^2 ~ ^2)52 — ^2) = 0- Hence, if we make 
the final choice of the base coordinates to be 

{el-ej){x-Y-ejz) ^ {el{x - ejz) - o\y) 
■' {ei-ei)z ' ^ x-Y-eiz 

the points tti, . . . , 7r4 will be on the line f + g = 0. Specifically, we have: 

TTl{dl+ Kl,-0l- Kl), TT3{9l+K3,-el-K3), 7r5(oO, 61^), TTy (S^ - , Oo) , 

TT2{el+ K2,~el~ K2), 7r4(0,0), TTe{oo,el), 7Ts{el + i,oo). 

Thus, we immediately get the identification between the parameters in the Riemann scheme of our Fuch- 
sian system and the parameters hi in the standard form of the Painleve equation: 

bi^el + Ki, b2^el + K2, 63 = ^? + Kg, 64 = 0, 65 = ^2. k^el b7 = e\-el, b^^-el-i. 

Then S ~ bi + ■ ■ ■ + bg, = —\ and the induced action of the standard Painleve dynamic on the Riemann 
scheme is 



a; = X ^ 1 X ^ 00 

el el Ki 

el el K2 

K3 



d-p(A(i)*) 
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whereas our elementary Schlesinger transformation acts as 

{}?} 






U2 


oo 


'1-1 


e\ + i 




el 




K2 











Thus, these two transformations correspond to the different translation directions in the symmetry root 
sub-lattice of the surface X and so are not equivalent. Indeed, we compute the action of ip^f of an elementary 
Schlesinger transformation { J ^ } on the classes "H/ , V-g, and £i to be 

Hf^2Hf + mg-£i-£2-£3-£i-'i£b- 2^8, 
ng^mf + bUg - 2£i - 2£2 - 2^3 - 2£i - - £& - 2£s, 

£i ^/ + 2'Hg — £2 — £z ^ £i ^ £b ^ 

£2 Hf + 2'Hg — f 1 — £3 — £4 — £5 — 

£3 1^ Hf + 2'Hg — f 1 — £2 — ^4 — ^^5 — ^81 

£4^nf + 2ng-£i-£2-£3-£5- £8, 
£5 '-^ £7, 

£6^2Hf + 2Hg -£i-£2-£3-£a- 2^5 - ^8, 
£T^2'Hf + 3Hg -£i-£2-£3~£a- 2^5 -£&- 2£i, 
£s '-^ Hg — £5, 

and compare with the standard dynamic ip^ given by (5.31-5.40) to see this explicitly: 

V-* : {ao,ai,a2,a3,a4,a5,ae) t-^ (ao, ai, 02, Q;3, "4, "5, "e) + (0, 0, 0, -1, 1, 1, -1)(5, 
LP* : (ao, ai,a2, ^3, a4,Q;5, ae) H> (ao, cti,a2, a3,a4,a5, "e) + (0,0,0, 1, 0, -1,0)5, 

where in both cases S — —Kx- 

In fact, wc can represent the standard Painleve dynamic as a composition of two different Schlesinger 
transformations, combined with automorphisms of our Fuchsian system, as follows. 

Let ai{j,k) be the map that exchanges the j-th and the fc-th eigenvectors of A^. If the eigenvalues ef 
and e^ are non-zero, this map is just a permutation on the decomposition space B x C. For this example we 
need ai(l,3), but it also defines a mapping on S x C, since bf e Ker(C|) and cf^ € Ker(Bi). Combining 
(Ti(l,3) with the scalar gauge transformation pi(— 0^), where pi{s) : A(x) i— (x — Mi)''A(a;), we get the action 
I]i(l, 3) = pi(— 0|[)o(7i(l, 3) on the spectral data (as encoded in the Riemann scheme) and the decomposition 
space: 



Ei(l,3) 




I --^l "^Z I T , 

<^i(i,3) n fli ^ Pii'Sl) 





Ei(l,3) : (hlhlcl\cf-M,hlcl\cl^) ^ (c^t x cf )S b^; " ; b^, b^; c 

^ ^ \ (c/ X )(b{ X b^j 



If „2t 
2 ' '-2 



where we had to take into account the normalization condition C|Bi = ©i. Then, looking at the action on 
the Riemann scheme, we see that 

d-P(A^^)*) = £1(1,3)0 {2 i}oEi(l,3)o {2 1}. 
Using Mathematica, we can verify directly that this composition dynamic satisfies equation (5.29). 

6. Conclusions 

To summarize, we gave an explicit expression for the generating function for elementary Schlesinger 
transformations in terms of the coordinates on the decomposition space for the coefficients of the underlying 
Fuchsian system. We also argued that this function should be thought of as a discrete Hamiltonian of 
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the dynamic generated by Schlesinger transformations, thus continuing to estabhsh a parallel between the 
theories of continuous and the discrete isomonodromic transformations. We used this discrete Hamiltonian 
to compute some explicit examples, including the derivation of the discrete Painleve equation of type d- 
P(^A^2^*'j in this framework, thus answering an earlier question of Sakai about the represent ability of this 
equation in terms of Schlesinger transformations. We also emphasized the role played by the geometry of 
rational surfaces corresponding to the equations in comparing different equations of the same type. 
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